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Abstract 

We argue that most commonly used models for nuclear scattering at ultra-relativistic 
energies do not treat energy conservation in a consistent fashion. Demanding theoretical 
consistency as a minimal requirement for a realistic model, we provide a solution for the 
above-mentioned problem, the so-called "Parton-Based Gribov-Regge Theory". 

In order to keep a clean picture, we do not consider secondary interactions. We provide a 
very transparent extrapolation of the physics of more elementary interactions towards nucleus- 
nucleus scattering, without considering any nuclear effects due to final state interactions. In 
this sense we consider our model a realistic and consistent approach to describe the initial 
stage of nuclear collisions. 
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1 Introduction 

The purpose of this paper is to provide the theoretical framework to treat hadron-hadron scattering 
and the initial stage of nucleus-nucleus collisions at ultra-relativistic energies, in particular 
with view to RHIC and LHC. The knowledge of these initial collisions is crucial for any theoretical 
treatment of parton thermalization and a possible parton-hadron phase transition, the detection 
of which being the ultimate aim of all the efforts of colliding heavy ions at very high energies. 

So what are the currently used models? Quite popular are semi-classical treatments of either 
partons or hadrons, in any case completely ignoring quantum mechanical interference. This is 
certainly quite unrealistic, and we do not want to discuss such options any further. 

There are as well considerable efforts to describe nuclear collisions via solving classical Yang- 
Mills equations, which allows to calculate inclusive parton distributions [Q. This approach is to 
some extent orthogonal to ours being based on the assumption of the perturbative nature of the 
triple Pomeron coupling The physical picture which stays behind the construction of Jl|] was 
outlined in |Q| and corresponds to perturbative interactions between individual parton cascades in 
hadrons, i.e. to fusion of partons of virtualities Q 2 > Qq, with Qq being a reasonable cutoff for 
QCD being applicable. Contrary to that, we believe that parton cascades interact with each other 
in the non-perturbative region of parton virtualities Q 2 < Qq and consider it as the interaction 
between soft Pomerons. In our scheme a relatively big value of the soft triple-Pomeron coupling 
should provide the necessary screening corrections which finally prevent the large increase of parton 
densities in the small x limit and restore the unitarity, thus leaving little room for higher twist 
effects in the perturbative part of the interaction. 

Another approach is the so-called Gribov-Regge theory (GRT) J|,||]. This is an effective field 
theory, which allows multiple interactions to happen "in parallel", with the phenomenological object 
called "Pomeron" representing an elementary interaction. Using the general rules of field theory, 
one may express cross sections in terms of a couple of parameters characterizing the Pomeron. 
Interference terms are crucial, they assure the unitarity of the theory. Here one observes an 
inconsistency: the fact that energy needs to be shared between many Pomerons in case of multiple 
scattering is well taken into account when calculating particle production (in particular in Monte 
Carlo applications), but energy conservation is not taken care of in cross section calculations. This 
is a serious problem and makes the whole approach inconsistent. Related to the above problem 
is the fact that different elementary interactions in case of multiple scattering are usually not 
treated equally, so the first interaction is usually considered to be quite different compared to the 
subsequent ones. 

Provided factorization works for nuclear collisions, one may employ the parton model, which 
allows to calculate inclusive cross sections as a convolution of an elementary cross section with 
parton distribution functions, with these distribution functions taken from deep inelastic scattering. 
In order to get exclusive parton level cross sections, some additional assumptions are needed, which 
follow quite closely the Gribov-Regge approach, encountering the same difficulties. 

As a solution of the above-mentioned problems, we present a new approach which we call 
"Parton-based Gribov-Regge Theory": we have a consistent treatment for calculating cross 
sections and particle production considering energy conservation in both cases; in addition, we 
introduce hard processes in a natural way, and, compared to the parton model, we can deal with 
exclusive cross sections without arbitrary assumptions. A single set of parameters is sufficient to fit 
many basic spectra in proton-proton and lepton-nucleon scattering, as well as in electron-positron 
annihilation (with the exception of one parameter which needs to be changed in order to optimize 
electron-positron transverse momentum spectra). 

The basic guideline of our approach is theoretical consistency. We cannot derive everything 
from first principles, but we use rigorously the language of field theory to make sure not to violate 
basic laws of physics, which is easily done in more phenomenological treatments (see discussion 
above) . 

There are still problems and open questions: there is clearly a problem with unitarity at very 
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high energies, which should be cured by considering screening corrections due to so-called triple- 
Pomeron interactions, which we do not treat rigorously at present but which is our next project. 

2 Problems 

Before presenting new theoretical ideas, we want to discuss the open problems in the parton model 
approach and in Gribov-Regge theory. 

Gribov-Regge Theory 

Gribov-Regge theory is by construction a multiple scattering theory. The elementary interactions 
are realized by complex objects called "Pomerons", who's precise nature is not known, and which 
are therefore simply parameterized: the elastic amplitude T corresponding to a single Pomeron 
exchange is given as 

T(s,t) ~is aa+alt 

with a couple of parameters to be determined by experiment. Even in hadron-hadron scattering, 
several of these Pomerons are exchanged in parallel, see fig. [j]. Using general rules of field theory 




Figure 1: Hadron-hadron scattering in GRT. The thick lines between the hadrons (incoming lines) represent a 
Pomeron each. The different Pomeron exchanges occur in parallel. 

(cutting rules), one obtains an expression for the inelastic cross section, 

4ct 2 = / dMl-cxp(-GM))}, (1) 

where the so-called eikonal G(s, b) (proportional to the Fourier transform of T(s, t)) represents one 
elementary interaction (a thick line in fig. [l]). One can generalize to nucleus-nucleus collisions, 
where corresponding formulas for cross sections may be derived. 

In order to calculate exclusive particle production, one needs to know how to share the energy 
between the individual elementary interactions in case of multiple scattering. We do not want to 
discuss the different recipes used to do the energy sharing (in particular in Monte Carlo applica- 
tions). The point is, whatever procedure is used, this is not taken into account in the calculation of 
cross sections discussed above. So, actually, one is using two different models for cross section cal- 
culations and for treating particle production. Taking energy conservation into account in exactly 
the same way will modify the cross section results considerably. 

This problem has first been discussed in ^,(^|. The authors claim that following from the non- 
planar structure of the corresponding diagrams, conserving energy and momentum in a consistent 
way is crucial, and therefore the incident energy has to be shared between the different elementary 
interactions, both real and virtual ones. 

Another very unpleasant and unsatisfactory feature of most "recipes" for particle production is 
the fact, that the second Pomeron and the subsequent ones are treated differently than the first 
one, although in the above-mentioned formula for the cross section all Pomerons are considered to 
be identical. 



3 A Solution: Parton-based Gribov-Regge Theory 
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The Parton Model 

The standard parton model approach to hadron-hadron or also nucleus-nucleus scattering amounts 
to presenting the partons of projectile and target by momentum distribution functions, /; u and 
fh 2 , and calculating inclusive cross sections for the production of parton jets with the squared 
transverse momentum p\ larger than some cutoff Qq as 

<t = E/ dpi J dx+J dx-fi 1 (x+,pl)fi 2 (x- )P l)^-(x+x- S )6( P l - Ql) , 

where dbijjdp\ is the elementary parton-parton cross section and i,j represent parton flavors. 

This simple factorization formula is the result of cancelations of complicated diagrams (AGK 
cancelations) and hides therefore the complicated multiple scattering structure of the reaction. 
The most obvious manifestation of such a structure is the fact that at high energies (y/s ^> 10 
GeV) the inclusive cross section in proton-(anti-)proton scattering exceeds the total one, so the 
average number N-^. of elementary interactions must be greater than one: 

T\Th 1 h 2 _ hih 2 I hifl2 v. i 
JV int — a mcl / a tot 1 ■ 

The usual solution is the so-called eikonalization, which amounts to re-introducing multiple scat- 
tering, based on the above formula for the inclusive cross section: 

^net 2 (s) = / rf 2 6{l-exp(-A(6)af n ^( S ))}=^a^' l Hs), (2) 

with 

a^{s) = J d%± ^exp (-A(b)a^( S )) (3) 

representing the cross section for n scatterings; A(b) being the proton-proton overlap function 
(the convolution of two proton profiles). In this way the multiple scattering is "recovered". The 
disadvantage is that this method does not provide any clue how to proceed for nucleus-nucleus (AB) 
collisions. One usually assumes the proton-proton cross section for each individual nucleon-nucleon 
pair of a AB system. We can demonstrate that this assumption is incorrect (see |Q). 

Another problem, in fact the same one as discussed earlier for the GRT, arises in the case of 
exclusive calculations (event generation) , since the above formulas do not provide any information 
on how to share the energy between many elementary interactions. The Pythia- method H amounts 
to generating the first elementary interaction according to the inclusive differential cross section, 
then taking the remaining energy for the second one and so on. In this way, the event generation 
will reproduce the theoretical inclusive spectrum for hadron-hadron interaction (by construction) , 
as it should be. The method is, however, very arbitrary, and - even more serious - we observe the 
same inconsistency as in the Gribov-Regge approach: energy conservation is not at all taken care 
of in the above formulas for cross section calculations. 



3 A Solution: Parton-based Gribov-Regge Theory 

In this paper, we present a new approach for hadronic interactions and for the initial stage of nuclear 
collisions, which is able to solve several of the above-mentioned problems. We provide a rigorous 
treatment of the multiple scattering aspect, such that questions as energy conservation are clearly 
determined by the rules of field theory, both for cross section and particle production calculations. 
In both (!) cases, energy is properly shared between the different interactions happening in parallel, 
see fig. H for proton-proton and fig. ||for proton-nucleus collisions (generalization to nucleus- nucleus 
is obvious). This is the most important and new aspect of our approach, which we consider to be 
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Figure 2: Graphical representation of a contribution to the elastic amplitude of proton-proton scattering. Here, 
energy conservation is taken into account: the energy of the incoming protons is shared among several "constituents" 
(shown by splitting the nucleon lines into several constituent lines), and so each Pomeron disposes of only a fraction 
of the total energy, such that the total energy is conserved. 




Figure 3: Graphical representation of a contribution to the elastic amplitude of proton-nucleus scattering, or more 
precisely a proton interacting with (for simplicity) two target nucleons, taking into account energy conservation. 
Here, the energy of the incoming proton is shared between all the constituents, which now provide the energy for 
interacting with two target nucleons. 

a first necessary step to construct a consistent model for high energy nuclear scattering. 

The elementary interactions, shown as the thick lines in the above figures, are in fact a sum 
of a soft, a hard, and a semi-hard contribution, providing a consistent treatment of soft and hard 
scattering. To some extend, our approach provides a link between the Gribov-Regge approach and 
the parton model, we call it "Parton-based Gribov-Regge Theory". 

4 Parton-Parton Scattering 

Let us first investigate parton-parton scattering, before constructing a multiple scattering theory 
for hadronic and nuclear scattering. 

We distinguish three types of elementary parton-parton scatterings, referred to as "soft", "hard" 
and "semi-hard", which we are going to discuss briefly in the following. The detailed derivations 
can be found in 

The Soft Contribution 

Let us first consider the pure non-perturbative contribution, where all virtual partons appearing in 
the internal structure of the diagram have restricted virtualities Q 2 < Qq, where Qq ~ 1 GeV 2 is a 
reasonable cutoff for perturbative QCD being applicable. Such soft non-perturbative dynamics is 
known to dominate hadron-hadron interactions at not too high energies. Lacking methods to cal- 
culate this contribution from first principles, it is simply parameterized and graphically represented 
as a 'blob', see fig. |j. It is traditionally assumed to correspond to multi-peripheral production of 
partons (and final hadrons) Q and is described by the phenomenological soft Pomeron exchange 
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Figure 4: The soft contribution. 

amplitude T so { t (s,t) |J. The corresponding profile function is expressed via the amplitude T ao ft as 
D so ft{s,b) = J d2 1i- exp(-iq±bj 21mT soft (,§, -q]_) 

27 P 2 art fsy-^- 1 f b 2 \ 

>&*/««.) u J exp r4A^ /so) ;' () 

with 

A ioft( Z ) = ni? part + "soft ln2, 

where s is the usual Mandelstam variable for parton-parton scattering. The parameters a so f t (0), 
ol soft are the intercept and the slope of the Pomeron trajectory, 7 pa rt and i?p art are the vertex value 
and the slope for the Pomeron-parton coupling, and sq ~ 1 GeV 2 is the characteristic hadronic 
mass scale. The external legs of the diagram of fig. |]are "partonic constituents", which are assumed 
to be quark-anti-quark pairs. 



The Hard Contribution 

Let us now consider the other extreme, when all the processes are perturbative, i.e. all internal in- 
termediate partons are characterized by large virtualities Q 2 > Qq- In that case, the corresponding 
hard parton-parton scattering amplitude can be calculated using the perturbative QCD techniques 
i fTol , Q| , and the intermediate states contributing to the absorptive part of the amplitude can 
be defined in the parton basis. In the leading logarithmic approximation of QCD, summing up 
terms where each (small) running QCD coupling constant a s (Q 2 ) appears together with a large 
logarithm ln(Q 2 /Aq CD ) (with Xqcd being the infrared QCD scale), and making use of the factor- 
ization hypothesis, one obtains the contribution of the corresponding cut diagram for t = q 2 = 
as the cut parton ladder cross section (T^ Td (s, Qq) as shown in fig. ]|, where all horizontal rungs 




Figure 5: The hard (or val-val) contribution. 

are the final (on-shell) partons and the virtualities of the virtual f-channel partons increase from 
1 Strictly speaking, one obtains the ladder representation for the process only using axial gauge. 
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the ends of the ladder towards the largest momentum transfer parton-parton process (indicated 
symbolically by the 'blob' in the middle of the ladder): 

= Kj2j dx+dx B dpl-^(x+x B s,pi) 

ml 

x £qcd(Qo, Ml, 4) E&aiQl M 2 , x B )9(M 2 - Q 2 ) , 

Here da'^ ln /dp 2 L is the differential 2 — > 2 parton scattering cross section, p\ is the parton transverse 
momentum in the hard process, m, I and Xg are correspondingly the types and the shares of the 
light cone momenta of the partons participating in the hard process, and Mp is the factorization 
scale for the process (we use Mp = Pj_/4). The 'evolution function' E^, T) (Qq, M F , z) represents 
the evolution of a parton cascade from the scale Qq to Mp, i.e. it gives the number density of 
partons of type m with the momentum share z at the virtuality scale Mp, resulted from the 
evolution of the initial parton j, taken at the virtuality scale Q 2 . The evolution function satisfies 
the usual DGLAP equation with the initial condition Eq^ u (Qq, Q\, z) = 8° m (5(1 — z). The 
factor K ~ 1.5 takes effectively into account higher order QCD corrections. 

In the following, we shall need to know the contribution of the uncut parton ladder r^ ard (s, t, Qq) 
with some momentum transfer q along the ladder (with t — q 2 ). The behavior of the corresponding 
amplitudes was studied in ||l3| in the leading logarithmic(l/x ) approximation of QCD. The precise 
form of the corresponding amplitude is not important for our application; we just use some of the 
results of namely that one can neglect the real part of this amplitude and that it is nearly 
independent on t, i.e. that the slope of the hard interaction i?j 2 lard is negligible small, i.e. compared 
to the soft Pomeron slope one has i?j; ard — 0. So we parameterize T^ ard (s , t , Qq) in the region of 
small t as Q 

T hard(*' Qa) = ig °har d ( J ! Qo) eX P (^hard *) (5) 

The corresponding profile function is obtained by calculating the Fourier transform Thard of 
Thard and dividing by the initial parton flux 2s, 

< rd (s,&) = ^2Imfij a fe rd (s,6), 

which gives 



DtA^-b) = ^rj d 2 g± exp(-iqj) 2JmT^ ld (s,- q l,Q 2 ) 

^-"'hard V ^hard / 

In fact, the above considerations are only correct for valence quarks, as discussed in detail in 
the next section. Therefore, we also talk about "valence- valence" contribution and we use D V al-val 
instead of Dhard: 

<i-valM)^ard(^)> 

so these are two names for one and the same object. 
The Semi-hard Contribution 

The discussion of the preceding section is not valid in case of sea quarks and gluons, since here 
the momentum share x\ of the "first" parton is typically very small, leading to an object with 
a large mass of the order Q\jx\ between the parton and the proton [jn|. Microscopically, such 



4 Parton-Parton Scattering 




Figure 6: The semi-hard "sea-sea" contribution: parton ladder plus "soft ends". 



'slow' partons with x\ <C 1 appear as a result of a long non-perturbative parton cascade, where 
each individual parton branching is characterized by a small momentum transfer squared Q 2 < Qq 
H [l(J. When calculating proton structure functions or high-pt jet production cross sections this 
non-perturbative contribution is usually included in parameterized initial parton momentum distri- 
butions at Q 2 = Qq. However, the description of inelastic hadronic interactions requires to treat it 
explicitly in order to account for secondary particles produced during such non-perturbative parton 
pre-evolution, and to describe correctly energy-momentum sharing between multiple elementary 
scatterings. As the underlying dynamics appears to be identical to the one of soft parton-parton 
scattering considered above, we treat this soft pre-evolution as the usual soft Pomeron emission, 
as discussed in detail in M. 

So for sea quarks and gluons, we consider so-called semi-hard interactions between parton 
constituents of initial hadrons, represented by a parton ladder with "soft ends", see fig. |[ As in 
the case of soft scattering, the external legs are quark-anti-quark pairs, connected to soft Pomerons. 
The outer partons of the ladder are on both sides sea quarks or gluons (therefore the index "sea- 
sea"). The central part is exactly the hard scattering considered in the preceding section. As 
discussed in length in [Q , the mathematical expression for the corresponding amplitude is given as 

jk 

with z being the momentum fraction of the external leg-partons of the parton ladder relative to 
the momenta of the initial (constituent) partons. The indices j and k refer to the flavor of these 
external ladder partons. The amplitudes Tj! oft are the soft Pomeron amplitudes discussed earlier, 
but with modified couplings, since the Pomerons are now connected to a parton ladder on one side. 
The arguments sq/z ± are the squared masses of the two soft Pomerons, z + z~s is the squared mass 
of the hard piece. 

Performing as usual the Fourier transform to the impact parameter representation and dividing 
by 2s, we obtain the profile function 

-^sea— sea ^) = 7T^ 21m T sca _ S ca{s- ^) > 

2s 

which may be written as 

Aea-sea (S, b) = £ / dz+ d Z ~ (z+) £*ft (*~) ^(z+zS, Qq) 

'^!'( .. . /' -I (7) 



4nX^ ft (l/(z+z-)) \ 4A£ ) ft (l/(z+z-)) 
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with the soft Pomeron slope A^ ft and the cross section o^ard being defined earlier. The functions 
E 3 sa{t (z^) representing the "soft ends" are defined as 



E j (z ± )=Jm.T :i (— 



t = 



■ > ("2) 

We neglected the small hard scattering slope i?£ ard compared to the Pomeron slope \ solt . We call 
Esoh also the " soft evolution", to indicate that we consider this as simply a continuation of the 
QCD evolution, however, in a region where perturbative techniques do not apply any more. As 
discussed in [Q, E 3 oit (z) has the meaning of the momentum distribution of parton j in the soft 
Pomeron. 

Consistency requires to also consider the mixed semi-hard contributions with a valence quark 
on one side and a non-valence participant (quark-anti-quark pair) on the other one, see fig. |. We 



^QCD 




Born 




Figure 7: Two "mixed" contributions. 



have 



,jk 
hard 



z s,q 



and 



D J val . sc AS,b) 



E 

k 



dz EL 



1 



•exp 



(8) 



where j is the flavor of the valence quark at the upper end of the ladder and k is the type of the 
parton on the lower ladder end. Again, we neglected the hard scattering slope R^ ard compared to 
the soft Pomeron slope. A contribution D 3 sca _ val (s, 6), corresponding to a valence quark participant 
from the target hadron, is given by the same expression, 

since eq. (||) stays unchanged under replacement z~ 
energy squared s for the parton-parton system. 



n val-sea ^) ' 



and only depends on the total cm. 



5 Hadron-Hadron Scattering 

To treat hadron-hadron scattering we use parton momentum Fock state expansion of hadron eigen- 
states jB| 

oo »\ k I k \ 

I^Efci/ l[dxirt(x 1 ,...x k )5[l- y £ f sj) a + ( Xl )---a+(x k )\0), 
k=i ■ J ° i=i \ j=i J 



5 Hadron-Hadron Scattering 



11 



where fk(%i, ■ ■ -Xk) ls the probability amplitude for the hadron h to consist of k constituent 
partons with the light cone momentum fractions x%,...,Xk and a + (x) is the creation operator for 
a parton with the fraction x, A general scattering process is described as a superposition of a 
number of pair-like scatterings between parton constituents of the projectile and target hadrons. 
Then hadron-hadron scattering amplitude is obtained as a convolution of individual parton-parton 
scattering amplitudes considered in the previous section and "inclusive" momentum distributions 
~\Fh^ { X U ■ ■ ■ x n) of n "participating" parton constituents involved in the scattering process(n > 1), 
with 

h P * n){Xl >--- Xn) = ^h n\(k-n)\ L II** l/*(*i,...**)| a *(l-£*. 

k — Th I — I - 1 \ J — 1 

We assume that Pj^hJxi, . . .x n ) can be represented in a factorized form as a product of the 
contributions Fp art (x/), depending on the momentum shares xi of the "participating" or "active" 
parton constituents, and on the function F^ emn ( 1 — Sj=i x j ) > representing the contribution of all 



spectators 



nucleon 




participants: 1 2 ... k 



Figure 8: Nucleon Fock state. 



"spectator" partons, sharing the remaining share 1 — J2j x j °f the initial light cone momentum (see 
% |): 

n I ii \ 

Pj i n \x 1 ,...x n )=Y\FX. t (x l ) F?ll-J2xA (9) 

The participating parton constituents are assumed to be quark-anti-quark pairs (not necessarily of 
identical flavors), such that the baryon numbers of the projectile and of the target are conserved. 
Then, as discussed in detail in the hadron-hadron amplitude may be written as 



no ^ J ti n [~l f 

iT hlh2 [s,t) = 87r 2 s^— / Y[dx+dxJ- f[ d 2 qi ± iT^ h ^ 

n=l U - J 1=1 Z=lL 07r Sl J 



+ - _ 2 



i - jytj F - m „^i - jyj <5 (2) (e&x - 9ij , (io) 

where the partonic amplitudes are defined as 

IP soft ' "'sea— sea ' val-val ' val— sea -'sea— val> 

with the individual contributions representing the "elementary partonic interactions plus external 
legs". The soft or semi-hard sea-sea contributions are given as 

r ^7«.-»a(^ + ^~.*.-9x) = T soit/sca - sca ( S ,-ql) F^ t {x+)F^ t {x-) 

x exp(-[R 2 hi +Rl 2 ]qi), (11) 
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the hard contribution is 

^-valC 3 ^ - ,^ 2 ) = / efe+^-p- [ dx-^^2T^ rd (x+ X - S ,q 2 ,Ql) 

Jo x vi j x v k 

x Fj$(x+,x+ - x+) F^(x~,x- - x-) exp(- [R 2 hi + R 2 h2 ] qfj , 
the mixed semi-hard "val-sea" contribution is given as 

Jo x v 

x F$g(x+,x + - x+) F^ It (x-) cxp(- [R 2 hi + R 2 h2 ] qf ± ) , 

and the contribution "sea-val" is finally obtained from "val-sea" by exchanging variables, 

T^^x+^s,?) = T% h J sea (x-,x+,s, q 2 ) . 

Here, we allow formally any number of valence type interactions (based on the fact that multiple 
valence type processes give negligible contribution). In the valence contributions, we have convo- 
lutions of hard parton scattering amplitudes T^ ard and valence quark distributions F^ art over the 
valence quark momentum share x^ rather than a simple product, since only the valence quarks are 
involved in the interactions, with the anti-quarks staying idle (the external legs carrying momenta 
x + and x~ are always quark-anti-quark pairs). 
The profile function 7 is as usual defined as 

lhih 2 {s,b) = — 2Imf/ tlfc3 (s,6), 
2s 

which may be evaluated using the AGK cutting rules with the result (assuming imaginary ampli- 
tudes) 

7^ 2 M) = E — / \\dx+dx-\{G% h ^x+,x-, Sl b) 

m=l J[J n=l £i=l 
1=0 ' J ° \=1 A=l 

^rcmn I ^ ^ ~~ ^ ^ ^\ J ^icmn I ^ ^ ^ ^ %\ J (-^) 

V a / V A / 

with a;P r °j/ tar s = 1 — J^xJ being the momentum fraction of the projectile/target remnant, and 
with a partonic profile function G given as 

G%*>(x+,x^a,b) = j-2lmf^(x+,x^,s,b), (13) 



see fig. g. This is a very important result, allowing to express the total profile function r )h 1 h 2 via 

'1TP 



the elementary profile functions G-, 91 2 



6 Nucleus-Nucleus Scattering 

We generalize the discussion of the last section in order to treat nucleus-nucleus scattering. In 
the Glauber-Gribov approach [jfj], 0], the nucleus-nucleus scattering amplitude is defined by the 



6 Nucleus-Nucleus Scattering 
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Figure 9: The hadronic profile function 7 expressed in terms of partonic profile functions G = G^ 



sum of contributions of diagrams, corresponding to multiple scattering processes between parton 
constituents of projectile and target nucleons. Nuclear form factors are supposed to be defined by 
the nuclear ground state wave functions. Assuming the nucleons to be uncorrelated, one can make 
the Fourier transform to obtain the amplitude in the impact parameter representation. Then, 
for given impact parameter bo between the nuclei, the only formal difference from the hadron- 
hadron case will be the averaging over nuclear ground states, which amounts to an integration over 
transverse nucleon coordinates bf and bf in the projectile and in the target respectively. We write 
this integration symbolically as 

/ dT AB := f JJ d 2 b? T A (b?) f[ d 2 bf T B (bf), (14) 
J J i=i j=i 

with A, B being the nuclear mass numbers and with the so-called nuclear thickness function T A (b) 
being defined as the integral over the nuclear density p A ( b) '■ 



T A (b) := J dz p A (b x ,b y ,z). (15) 

It is convenient to use the transverse distance bk between the two nucleons from the k-th nucleon- 
nucleon pair, i.e. 

bk = b Q + b* [k) - 6f (fc) 

where the functions 7r(fc) and r(fc) refer to the projectile and the target nucleons participating 
in the k th interaction (pair k). In order to simplify the notation, we define a vector b whose 
components are the overall impact parameter b as well as the transverse distances 61, ...,b A B of 
the nucleon pairs, 

b = {60,61, -.,b AB }- 

Then the nucleus-nucleus interaction cross section can be obtained applying the cutting procedure 
to elastic scattering diagram and written in the form 

0to3(«) = J d2b o J dT ABlAB (s,b), (16) 

where the so-called nuclear profile function j AB represents an interaction for given transverse 
coordinates of the nucleons. 

The calculation of the profile function ^ AB as the sum over all cut diagrams of the type shown 
in fig. |l(] does not differ from the hadron-hadron case and follows the rules formulated in the 
preceding section: 

• For a remnant carrying the light cone momentum fraction x (x + in case of projectile, or x~ 
in case of target), one has a factor F rcmn (x). 
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Figure 10: Example for a cut multiple scattering diagram, with cut (dashed lines) and uncut (full lines) elementary 
diagrams. This diagram can be translated directly into a formula for the inelastic cross section (see text). 

• For each cut elementary diagram (real elementary interaction = dashed vertical line) at- 
tached to two participants with light cone momentum fractions x + and x~ , one has a factor 
G(x + , x~ , s, b). Apart from x + and x~ , G is also a function of the total squared energy s 
and of the relative transverse distance b between the two corresponding nucleons (we use G 
as an abbreviation for G^ for nucleon-nucleon scattering). 

• For each uncut elementary diagram (virtual emissions = full vertical line) attached to two par- 
ticipants with light cone momentum fractions x + and x~ , one has a factor — G(x + , x~ , s, 6), 
with the same G as used for the cut diagrams. 

• Finally one sums over all possible numbers of cut and uncut Pomerons and integrates over 
the light cone momentum fractions. 

So we find 




with 

*r j = i- E x L, 

7r(k)—i 

*r = 1 e *w 

r(k)=] 

The summation indices mj~ refer to the number of cut elementary diagrams and Ik to number of 
uncut elementary diagrams, related to nucleon pair k. For each possible pair k (we have altogether 
AB pairs), we allow any number of cut and uncut diagrams. The integration variables x k refer to 
the /i th elementary interaction of the fc th pair for the cut elementary diagrams, the variables x^ x 
refer to the corresponding uncut elementary diagrams. The arguments of the remnant functions 
-Frcmn are the remnant light cone momentum fractions, i.e. unity minus the momentum fractions 



6 Nucleus-Nucleus Scattering 
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of all the corresponding elementary contributions (cut and uncut ones). We also introduce the 
variables xf IOJ and x^ 6 , denned as unity minus the momentum fractions of all the corresponding 
cut contributions, in order to integrate over the uncut ones (see below). 

The expression for jab sums up all possible numbers of cut Pomerons nik with one exception 
due to the factor (1 — Sos mk ): one does not consider the case of all m^'s being zero, which 
corresponds to "no interaction" and therefore does not contribute to the inelastic cross section. We 
may therefore define a quantity jab, representing "no interaction", by taking the expression for 
jab with (1 - S Sm k ) replaced by (<5o£m fc ): 

p AB f l k 1 AB r l h \ 

1AB ( S ,b) = / n n^A^A n n n-^^wA) 

h l AB k=l Ia=1 J fc=l I fe ' A=l J 



A / \ B 



X 



U F+ 1 E '/.a IP i E (is) 

»=1 V 7r(fc)=i / J=l \ r(*)=j 



One now may consider the sum of "interaction" and "no interaction", and one obtains easily 

j AB (s,b)+j AB (s,b) = l. (19) 

Based on this important result, we consider jab to be the probability to have an interaction and 
correspondingly jab to be the probability of no interaction, for fixed energy, impact parameter 
and nuclear configuration, specified by the transverse distances bk between nucleons, and we refer 
to eq. ( |l9|) as "unitarity relation". But we want to go even further and use an expansion of jab in 
order to obtain probability distributions for individual processes, which then serves as a basis for 
the calculations of exclusive quantities. 

The expansion of jab in terms of cut and uncut Pomerons as given above represents a sum of 
a large number of positive and negative terms, including all kinds of interferences, which excludes 
any probabilistic interpretation. We have therefore to perform summations of interference contri- 
butions - sum over any number of virtual elementary scatterings (uncut Pomerons) - for given 
non- interfering classes of diagrams with given numbers of real scatterings (cut Pomerons) . Let 
us write the formulas explicitly. We have 



lAB{s,b) 



. AB f m k ~| 

E-Ed-^J / II \U d <^A 

mi rn AB k=\ 1/1=1 ) 

AB ( 1 m k \ 

x II \—i II G(x + k ^x-^s,b k ) ®ab (x^,x^,s,b) , (20) 

k=l K " n=l ) 

where the function $ representing the sum over virtual emissions (uncut Pomerons) is given by 
the following expression 



„ ab r i k \ AB ( l k ~| 

$ AB (x^,x^, s ,b) = E -E / n n n n- G «A^A,*A) 

h l AB J k=l Ia=1 J k=l { k ' \=1 ) 

x nF remn (xr j - E -u)f[F^[ x ^ s - ^,a)- ( 21 ) 

i = l \ 7T(fc)=i / j=l \ T(k)=j I 



This summation has to be carried out, before we may use the expansion of jab to obtain probability 
distributions. This is far from trivial, the necessary methods are described in [Q. To make the 
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notation more compact, we define matrices X + and X , as well as a vector to, via 



X + = 

x- = 



R} . 

{mk}, 



which leads to 



lAB {s.b) = J2( l ~ S o™) J dX + dX~n ( ^\m,X+,X-), 



lAsis^b) = n^P (0,0,0), 



with 



proj targ 



1 ^ 1 



,,b) 



AB ( m k -\ 

rfjg\m,X+,X-) = ]l\—I[G(x+^,8,b k ) \ <t> AB ( 

k=l { M =l J 

This allows to rewrite the unitarity relation eq. ( |l9| ) in the following form, 

[ dX+dX-Vl { ^\m,X + ,X-) = 1. 

This equation is of fundamental importance, because it allows us to interpret fl( s ' b >(m,X + , X~) 
as probability density of having an interaction configuration characterized by to, with the light 
cone momentum fractions of the Pomerons being given by X + and X~ . 

7 Virtual Emissions and Markov Chain Techniques 

What did we achieve so far? We have formulated a well denned model, introduced by using the 
language of field theory, solving in this way the severe consistency problems of the most popular 
current approaches. To proceed further, one needs to solve two fundamental problems: 

• the sum $ab over virtual emissions has to be performed, 

• tools have to be developed to deal with the multidimensional probability distribution fi)^ , 

both being very difficult tasks. Introducing new numerical techniques, we were able to solve both 
problems, as discussed in very detail in Qj]. 

Calculating the sum over virtual emissions ($ab) is done by parameterizing the functions G 
as analytical functions and performing analytical calculations. By studying the properties of &ab, 
we find that at very high energies the theory is no longer unitary without taking into account 
screening corrections due to triple Pomeron interactions. In this sense, we consider our work as a 
first step to construct a consistent model for high energy nuclear scattering, but there is still work 
to be done. 

Concerning the multidimensional probability distribution fi^£^ (to, X + , X~), we employ meth- 
ods well known in statistical physics (Markov chain techniques) . So finally, we are able to calculate 
the probability distribution rt^\m, X + , X~), and are able to generate (in a Monte Carlo fashion) 
configurations (m,X + ,X~) according to this probability distribution. 



8 Summary 
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8 Summary 

What are finally the principal features of our basic results, summarized in eqs. ([l(], ||(], |2l|)? 
Contrary to the traditional treatment (Gribov-Regge approach or parton model), all individual 
elementary contributions G depend explicitly on the light-cone momenta of the elementary in- 
teractions, with the total energy-momentum being precisely conserved. Another very important 
feature is the explicit dependence of the screening contribution <&ab (the contribution of virtual 
emissions) on the remnant momenta. The direct consequence of properly taking into account 
energy-momentum conservation in the multiple scattering process is the validity of the so-called 
AGK-cancelations in hadron-hadron and nucleus-nucleus collisions in the entire kinematical region. 

The formulas (|jl| ^O], |2lj) allow to develop a consistent scheme to simulate high energy nucleus- 
nucleus interactions. The corresponding Monte Carlo procedure is exactly based on the cross 
section formulas so that the entire model is fully self-consistent. 
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